Field-induced incommensurate order for the quasi-one-dimensional XXZ model in a 

magnetic field 
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We investigate phase transitions of the quasi-one-dimensional S = 1/2 XXZ model in a magnetic 
field, using the bosonization combined with the mean-field treatment of the inter-chain interaction. 
We then find that the field induced incommensurate order is certainly realized in the low field region, 
while the transverse staggered order appears in the high field region. On the basis of the result, we 
discuss the field- induced phase transition recently observed for BaCo2V2 08. 
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I. INTRODUCTION 

Field-induced phase transition in the quantum spin 
systems has been providing interesting physics such as 
magnon Bose Einstein condensationi. Recently an ex- 
otic field induced phase transition was observed for 
BaCo2V2 0g'^, which can be regarded as the quasi-one- 
dimensional(lD) S — 1/2 XXZ antiferromagnt hav- 
ing the Ising-hke anisotropy A ~ 2; the magnetization 
and electron spin resonance (ESR) measurements above 
1.8K show that BaCo2V2 08 is basically described by the 
Bethe-ansatz-based theoretical analysisi^ However, the 
specific heat measurements up to 12T below 1.8K has 
revealed that the weak 3D couplings possibly trigger the 
exotic incommensurate(IC) order in the low-field region."* 
A peculiar point on this phase is that the ordering is dif- 
ferent from the Neel type at the zero magnetic field and 
the shape of the phase boundary in the H-T plane is 
quite different from the usual field induced order in the 
coupled Haldane system^. This suggests that the Ising- 
like anisotropy in the quasi- ID system plays an essential 
role in the field induced IC order phase, behind which 
there is substantially important physics. 

The ID XXZ antiferromagnet is an exactly solved 
model playing the essential role to understand the critical 
quantum fluctuation and strong correlation effects. 
Although the Ising-like anisotropy favors the z directed 
Neel(z-Neel) order at the zero magnetization, the mag- 
netic field beyond the critical field He recovers the critical 
quantum fluctuation(see FiglT]) and then the system is de- 
scribed by the Tomonaga-Luttinger(TL) liquidji which is 
characterized by the power law decay of the correlation 
f mictions: 
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where rj is the TL exponent, Mz is the uniform magneti- 
zation due to a magnetic field H , and the corresponding 
Fermi wave number is kp = '"'(1/2 — Mz). The nonuni- 
versal coefficients A and B were evaluated in Refi^. For 
the isotropic Heisenberg model, < 1 is always satisfied 



and thus the transverse fluctuation of ^ is dominant. 
For the Ising-like case, however, 77 > 1 appears in the 
low field region, where the longitudinal IC fluctuation 
becomes dominant. In the actual quasi-lD compound, 
there is indispensable inter-chain interaction, which may 
bring the finite temperature phase transition accompa- 
nying the field-dependent IC order. 

In this paper, we investigate the field induced IC order 
for the coupled XXZ chains, using the bosonization com- 
bined with the mean-field treatment of the inter-chain 
interaction*''. In particular, we make quantitative analy- 
sis of the transition temperatures, taking account of the 
non-universal coefficients in ^ and Wc then find 
that the IC order is certainly realized in the low field 
region, in contrast with Refii where the possibility of 
the IC order was not taken into account, while the trans- 
verse staggered order appears in the high-field region. 
Moreover, we show that the present theory successfully 
explains the field dependence of the experimentally ob- 
served transition temperature. We also determine the 
inter-chain coupling of the BaCo2V2 08 as 0.09K. 

This paper is organized as follows. In Sec. II, we ex- 
plain the model and the mean-field theory for the inter- 
chain interaction on the basis of bosonization. In Sec. Ill, 
magnetic-field dependences of the transition tempera- 
tures are presented for the IC order and the transverse 
staggered order. Then the IC order of BaCo2V208 is dis- 
cussed in detail. In Sec. IV, we summarize conclusions 
and discuss related topics. 



II. MODEL AND FORMULATION 

The relevant model we consider here is the weakly cou- 
pled S — 1/2 XXZ chains on the simple cubic lattice, 
whose Hamiltonian is given by 

W = '^^[J[Sn,j ■ Sn+l,j)A 
n.j 

+ J2 J'{S^,rSn,r)A^HY,Sl,, (3) 

where [S ■ S)a = S^'S'' + S'^Sy + A^^S*^ is the deformed 
inner product. The index n runs along the chain direc- 



2 




FIG. 1: (online color) The exact magnetization curve M, crit- 
ical exponent r] and spin wave velocity v for the XXZ chain 
of A = 2. 



tion, j labels the inter-chain directions, and 0, j') de- 
notes the nearest neighbor pair of the chains. Then J is 
the exchange coupling along the chain direction and the 
inter-chain interaction is controlled by J'{<^ J). Note 
that we set the lattice space to be unity. 

Let us discuss the order-disorder transition in _ff > He- 
Since the simple cubic lattice is considered, we can set 
up the "sub-chain" mean fields as Sn,j = M„ + SSnj, 
where Mn is a classical vector field. We assume the 
IC oscillation of the magnetization of the z component 
around the uniform(average) magnetization, 



[Mz ± rriic cos(2fci?n)]e2 



(4) 



where rriic is the amplitude of the IC fluctuation around 
the average magnetization and the sign depends on 
the sub chain. Also for the transverse staggered fluctua- 
tion, we assume 



(5) 



where is the amplitude of the staggered magnetiza- 
tion in the transverse direction. Then, the mean-field 
treatment of the inter-chain coupling yields the mean- 
field Hamiltonian per chain as 



ic/s 



where 
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is the ID XXZ chain in an effective magnetic field H 
and 7i^c/s perturbation originating from the mean 

fields. For the IC order, we have H = H — zAJ'AL and 



. cos{2kpn)S^ + const. 



(8) 



where hie = zAJ'm. Note that the coordination number 
is z = 4 for 3D. These effective fields should be deter- 
mined self-consistently. Also for the transverse staggered 



fluctuation, we have H = H ~ zA J'Af, and 



H, = -/i.2^(-)"5^^+ const 



(9) 



with hg = zj'uis 

In order to treat the IC nature in the mean-field Hamil- 
tonian, it is useful to employ the effective model in 
the continuum limit. This is achieved by the standard 
bosonization scheme. Assuming H — He ^ hie or h^, 
we can write the XXZ chain in the magnetic field as 



Hid ^ 2 / ^^[(^-'^) 



(10) 



where v is the spin wave velocity and the compactification 
radius R is defined by the boundary condition (f)(x + L) = 
4>{x) + ttRN. The equal time commutator of the bosonic 
fields is defined as [(j){x),9{y)] ~ iQ{x—y), where Q{x—y) 
is the step function. For the XXZ model, the radius varies 
R = I/a/tF —> l/\/47r as the magnetic field increasing 
from He to the saturation(free fermion) limit. The TL 
exponent is given by 77 = 2TrR^. 

The boson representation of the spin operators is given 
by the formula 
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where a and b are nonuniversal coefficients depending on 
A and H. Using ([TT]) and ([12]), the amplitude of the 
equal time correlation function is given by Ai — a^/2, 
Bq ~ b'^/2 with the rcgularization 



dk- 
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(13) 



where a is a cutoff parameter. Although the analytical 
expression of these coefficients in the magnetic field are 
still unknown, the numerical value is available in Ref.— , 
which play a crucial role to semi-quantitative evaluation 
of the transition temperature, as will be seen later. Sub- 
stituting ([TT|) and ([T2I) into ^ and (O, we obtain the 
boson field representation of the perturbations. 
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where we have omitted the 2kF and Akp oscillating terms. 

Let us consider the effects of Ti-ic/s on the Hamiltonian 
PH)) . The gap generated by the staggered field was 
analyzed in Refji^ to be AEs hV^^~^^ ■ Also for the 
operator of cos(0/i?), the similar standard renormaliza- 
tion group argument leads AEie ^hc^^"^^^^ ■ Both of 
the operators are always relevant between He and sat- 
uration field. In the lower field region, however, 77 > 1 
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and thus 2/(4 - 1/ry) < 2/(4 - 77); The IC field is more 
relevant in the low field region, while the transverse stag- 
gered field is more relevant in the high field region. The 
border is just 77 = 1 namely, an effective SU(2) point. 
For the case of A = 2, it corresponds to H/ J ~ 1.5 as 
in Figdl This analysis of the gap is consistent with the 
naive expectation of the critical exponent 77, which leads 
the IC order in the low field region. For the quantita- 
tive analysis, however, the coefficients of the gaps become 
important. It should be also noted that the assumption 
H — He 3> hie, hg is not valid in the vicinity of the lower 
critical field. 

Now we proceed finite temperature behaviors. In the 
framework of the mean-field theory, hie or hs should be 
determined self-consistently. Write the IC magnetization 
of ^ aXT , H and hie as fie{T, H, hs), and the staggered 
magnetization at T, H and hg as fs{T, H,hs). The self- 
consistent equations are written as m.^c — fie(T, H,hie) 
and rUg = fs(T, H,hs), combined respectively with hie — 
zAJ'rriie and hg = zJ'mg. Taking hie, hg limits, we 
can determine the transition temperatures 
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where Xic = dftc/dhie\hi,=o and Xs = dfg/dhg\h,=o. 
According to the linear response theory, the 
dynamical susceptibility can be represented 
through the correlation function: Xaf3{Q,'-^',T) = 
-i Y.n I die''^*~'9"e(t)([5"(n, t),S^{0, 0)])t, where 
(• • •)t denotes the average at a temperature T. For the 
TL Hamiltonian (jlOp . this dynamical susceptibility was 
actually calculated in Ref.™—* . For the estimation of 
the transition temperature, the susceptibility at the soft 
mode is essential: for the IC order, ui = and q = 2fci?, 
and for the staggered order, w = and q = tt. Taking 
account of p^ . the explicit form of the leading term 
becomes 



Xic = Xzz{2kF,0;T) = 



Tt\ /2^r\'/" ' 1 1 ^ 
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where B{x, y) is the Euler's beta function. Note that the 
cutoff parameter in p7p and ^TE\i formally corresponds 
to a = 1, due to the regularization (|13p . The magnetic- 
field dependence is implicitly included in v, 77, Ai and 
Bq. Substituting the above susceptibilities into PB)) . we 
obtain the final result of the transition temperatures 



Ti-) = ^(zAfA^^Bi'l-^fY^W) 
" 27r V V Ari' 2r]' J ^ ' 



(20) 



In the above expression, r] and v can be exactly calculated 
by solving the Bethe ansatz integral equation as in Fig. 
[1] In addition, the non universal amplitudes Ai and Bq 
can be obtained by using density matrix renormalization 
group and the bosonization expression of the correlation 
function for the open boundary system-. We can thus 
calculate Tc semi-quantitatively without any additional 
parameter. Here it should be noted that the previous es- 

(s) 

timation for Tc was based on the correlation amplitude 
at the zero magnetic field for A < 1—. 



III. RESULTS 
A. phase diagram 

On the basis of p9p and ((20|) . we calculate the mag- 
netic field dependences of the transition temperatures of 
J' = 0.01 for various A. The correlation amplitudes 
are extracted from the correlation functions obtained via 
density matrix renormalization group, as mentioned in 
the previous section. In figure 2, we show the result- 
ing phase diagrams in the T — H plane, where the solid 
and broken lines respectively indicate (fTQ]) and ([20|l . and 
the magnetic field is normalized by the saturation field 
Hg . The curve corresponding to the higher Te is realized 
as an actual order-disorder transition. In the following, 
we concentrate on the order-disorder transitions between 
the critical field He and the saturation field Hg. Thus 
the z-Neel phase of M — below He is not shown here 
explicitly. In addition, note that He for A = 1.05 and 
1.5 is in vicinity oi H ~ in the scale of FigO 

In Fig[21 we can see that the IC order certainly occurs 
above the critical field He- For A — 1.05, the IC order 
appears in the vicinity oi H = He — 0- As A increases, 
the transverse staggered order is suppressed, while the IC 
order develops rapidly and the corresponding range of H 
extends to the higher field region. An important feature 

(ic) 

of the IC order is that the field-dependence of Tc il- 
lustrates a characteristic curve: tI^""^ has the maximum 
near He, and it decreases rapidly as H increases. Such 
shape of the phase boundary should be contrasted to the 
semicircle-like boundary for the field-induced staggered- 
order in the coupled Haldane system. As further increas- 
ing H, the curves for T^*'^'' and ri''"' intersect at a certain 
magnetic field, which is denoted as i/* henceforth, and 
the transverse staggered order appears for Hg > H > H* . 
We can also see that, as A becomes large, H* shifts to 
the higher field side. 

The behaviors above are basically consistent with the 
argument based on the TL exponents of the XXZ chain, 
since the region of 77 > 1 appears above the critical field 
He and it extends rapidly to the higher field side, as A 
is increased. However, it should be noted that H* does 
not coincide with the effective SU(2) point ?y = 1. This is 
because 7; = 1 is achieved in the effective field theory level 
and thus Ai ^ Bq is permitted even at 77 = 1, in contrast 
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to the isotropic Heisenberg chain at the zero field having 
the SU(2) symmetry in the spin operator level. Since the 
correlation amplitude Ai has a larger value than i?o(e.g. 
see FIG. 2 in Ref.^), ri**^-* is relatively enhanced than 

(s) 

Tc , implying that H* slightly shifts to the higher field 
side than the effective SU(2) point. In this sense, the 
precise amplitudes are essential in the inter-chain mean- 
field theory. In addition, we can see that A in (fTO)) is also 
a source of such an enhancement of the IC order. 

We next discuss the inter-chain-coupling dependence 
of transition temperatures. According to eqs. (fTO]) and 

(PO)) . the precise J'-dependences are given by Tc^'^^ oc 
j/,,/(2,,-i) j,{s) ^ j/i/(2~,,)^ gQ ^^^^ ^j^g g(,j^lg ^Yie 

transition temperature naturally becomes large, as J' is 
increased. In particular, we can see that T^'"^ IS more eas- 
ily lifted toward the higher field region where rj < 1, since, 
as mentioned above, the IC order is basically enhanced by 
the correlation amplitude Ai and the anisotropy A. How- 
ever, we should remark that such enhancement of Ti*"^-* in 
the inter-chain mean-field theory does not always lead to 
a clear observation of the IC order for a larger J'; we need 
to pay special attention to stability of the IC order. Let 
us recall that the spin flop transition occurs for the case of 
the spatially isotropic exchange coupling(J = J'), where 
the magnetization directly jumps from z-Neel phase of 
M = to the transverse staggered ordered state. This 
suggests that the IC order becomes thermodynamically 
unstable beyond a certain critical J' so that it is em- 
bedded in the magnetization jump. Unfortunately, the 
critical J' can not determined within the frame work of 
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FIG. 2: (online color)The transition temperatures for A = 
1.05, 1.5, 2.0 and 3.0 witli J' = 0.01. Tlie intra-chain cou- 
pling J is set to be unity. The solid line and the broken 



line respectively represent ri*'^' and Tc''' ■ The vertical dotted 
lines for A = 2.0 and 3.0 indicate the critical field He- He 
for A = 1.05 and 1.5 is not shown here, since it is located in 
vicinity of H = 0. The horizontal axis is normalized by the 
saturation field Hs. 



the mean-field theory for the inter-chain coupling, since 
analytical calculation of the free energy is still a difficult 
task. In the next subsection, nevertheless, we shall show 
that the stable IC order actually occurs in the experi- 
mental situation of BaCo2V208. 



B. comparison with experiment 

Let us discuss the field dependence of the transi- 
tion temperature of BaCo2V208. The basic parame- 
ters were determined by the magnetization and ESR 
measurements^. The exchange coupling in the chain di- 
rection is given by J ~ 65K and the precise anisotropy 
parameter is A ~ 2.17. The critical field is He — 3.9T 
and the saturation field is about 23T with g — 6.2. In ad- 
dition, we have actually calculated the correlation ampli- 
tudes of A = 2.17. In Figure[31 the field dependence of the 
transition temperature for J' ~ 0.09K(J7J = 0.00138) 
is illustrated together with the experimental data, where 
the solid and broken lines indicate and (PO]) respec- 
tively. 
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FIG. 3; (online color)The transition temperatures for the 
IC order and the transverse staggered order. The solid line 
and the broken line respectively represent T^'^^ and Ti"' for 
A = 2.17 and J' / J = 0.00138. The vertical dotted line means 
He = 3.9T. The solid circles indicate the experimentally ob- 
served transition temperature'* and the solid triangles mean 
the phase boundary between the z-Neel order and the field 
induced IC order. 

In the left side of corresponding to the solid tri- 
angles, the z-Neel order occurs, at which the uniform 
magnetization is zero. Note that the transition temper- 
ature to the z-Neel phase at iJ = is about 5.4K, which 
is much higher than Tc^'^\ Above the z-Neel order is 
destroyed by the magnetic field and then we come into 
the targeted region of the present theory. The solid cir- 
cles indicate the experimental transition temperature up 
to 12T. We can see that the theoretical curve excel- 
lently reproduces the experimental results, implying the 
inter-chain mean-field theory is basically correct for 3D. 
A remarkable point is that the shape of the experimental 
phase boundary is consistent with the theoretical curve 
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of the IC order; As H increases above He, Tc decreases 
rapidly from T ~ 1.7K down to 0.4K. The above facts 
support that the IC order driven by the one dimension- 
ahty can be thcrmodynamically stabihzcd in the exper- 
imental situation. Another interesting point is that, as 

(ic) 

further increasing the theoretical curves for n ' and 

is) 

Tc intersects at H* ~ 15. IT, which predicts that the 
transverse staggered order appears above H ~ 15. IT. In 
order to verify the theory, a specific heat measurement 
in the higher field is is highly desirable. However, the 

(s) 

value of Tc is relatively low and thus the experimental 
observation in the competing region may be subtle. 

IV. SUMMARY AND DISCUSSIONS 

We have discussed the field induced IC order on the 
basis of the bosonization combined with the mean-field 
theory for the inter-chain interaction. In particular, the 
numerically exact correlation amplitudes plays the cru- 
cial role to explain the shape of the experimental phase 
boundary. In order to investigate the IC order beyond the 
mean-field level, we have also performed quantum mon- 
tecarlo (QMC) simulations based on the directed loop 
algorithm. Then, we have confirmed that the IC order 
actually occurs for a J'/ J = 0.1—. We can therefore con- 
clude that the field-induced IC order is certainly realized 
in the actual system and the inter-chain mean-field treat- 
ment captures the essential nature of it. The inter-chain 
coupling of BaCo2V208 estimated within the mean-field 
theory is J' ~ 0.09K. 

From theoretical point of view, the phase transition for 
the 3D classical spin model with the easy-axis anisotropy 
was intensively studied in 70s, in the context of the spin 
flop transitioniii The low magnetized state is unstable 
in the 3D isotropic lattice system and the magnetization 
jumps directly from the z-Neel state into the spin flopped 
state. The spin flop transition also occurs for the 2D 
Ising-like XXZ model on the isotropic square lattice at 
the zero temperature^^ The present result implies that, 
as the ID fluctuation is enhanced, the IC order — spin 
version of the charge density wave(CDW) — emerges in 
the phase diagram. Of course, BaCo2V2 08 is insulat- 



ing, and thus the mechanism is attributed to the nesting 
of "spin" itself. In this sense, the present IC order is 
very similar to that in the spin-Peierls systemi^^ How- 
ever, the driving mechanism is the inter-chain spin-spin 
interaction itself rather than a spin-phonon coupling in 
the spin-Peierls case. Since the inter-chain interaction fa- 
vors the transverse staggered order as well, the spin flop 
transition may be induced with a certain finite inter-chain 
coupling, implying that the thermodynamic stability of 
the IC order is a non-trivial question. The present result 
demonstrate that the IC order based on the ID mech- 
anism is certainly stabilized in the actual experimental 
situation. For the quasi- ID spin model, the Fermi wave 
number kp can be easily controlled by the magnetic field, 
in contrast with the CDW in the metallic system. A fur- 
ther experimental study, particularly neutron scattering 
experiment in the magnetic field, is highly interesting. 
In addition, the connection to the spin flop transition in 
the high field region is also theoretically important prob- 
lem, although the experiment for BaCo2V208 suggests 
a weak first order transition at He accompanying the 
spin-lattice coupling which may cooperatively stabilize 
the incommensurate order. 

Finally we remark that our theory is valid not only 
for the similar quasi- ID systems with the easy axis 
anisotropy, but also for a class of the frustrating systems. 
In fact, the frustrating systems is often mapped into an 
effective XXZ model, for which the IC order is actually 
pointed outi^i^i Such enhancement of the IC fluctua- 
tion is also reported for an anisotropic S = 1 chain22i. 
We hope that the rich physics associated with the spin 
anisotropy and quantum fluctuation can be developed by 
further theoretical and experimental researches. 
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